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Abstract. Wc study the low-lying spectrum of the Dirichlet Laplace operator on a ran- 
domly wiggled strip. More precisely, our results are formulated in terms of the eigenvalues 
of finite segment approximations of the infinite waveguide. Under appropriate weak- 
disorder assumptions we obtain deterministic and probabilistic bounds on the position of 
the lowest eigenvalue. A Combes-Thomas argument allows us to obtain a so-called 'initial 
length scale decay estimates' at they are employed in the proof of spectral localization 
using the multiscale analysis method. 

I. Introduction 

The propagation of waves in disordered media can be modeled using differential equations 
governed by a random Hamiltonian. The most important questions in this context concern 
the long time behavior of waves which oftentimes allows to conclude results concerning the 
transport properties of the material described by the random operator. A particularly well 
studied class of operators is the one which arises in the quantum mechanical description of 
disordered solids. To this class belong various types of random Schrodinger operators, e.g. 
the Anderson model on £ 2 (Z d ), or a Laplacian with Poisson distributed repulsive impurity 
potentials on L 2 (R d ). 

For most operators it is possible to relate propagation properties to spectral features, 
e.g. by the use of the RAGE or Ruelle theorem. From this point of view it is justified 
to study first the measure theoretic spectral types which arise in a certain random model, 
and then relate them to the transport properties of the considered material. This leads to 
the study of a plethora of spectral features, some of which are specific to the situation that 
we are dealing with a random operator, i.e. with an infinite family of individual operators. 
Let us mention some of these properties: the characteristics of the spectrum as a subset of 
the real line, its band structure, location of the spectral innmum (i.e. the innmum of the 
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bottom of the spectrum when varied over all members of the family), the density of the 
spectrum in various energy regions, decay properties of the Green's function, etc. 

To be able to point out the interesting contribution of the paper at hand, we shall assume 
in the further discussion that the randomness enters the Hamiltionan via a countable family 
of random variables. Already when considering the very basic features of the spectrum, one 
sees that it makes a great difference whether the operator (more precisely: the associated 
quadratic form) depends in a monotone or non-monotone way on the random variables 
entering the model. In the case of monotone dependence the spectral minimum of the 
operator family obviously corresponds to the configuration where all random variables 
are set to one of the extremal values. Similarly, in the monotone situation the band 
structure of the spectrum can be analyzed using rather basic sandwiching arguments, 
see e.g. |KSS98j . It is consistent with these elementary examples of the advantages of 
monotonicity that there is a rather good understanding of typical energy /disorder regimes 
where monotone models exhibit localization of waves, see the monographs and survey 
articles [StoOTT IKM071 IVes07l lKir08] . 

If the dependence of the quadratic form on the random variables is not monotone, already 
the identification of the spectral minimum is not obvious and sometimes a highly non-trivial 
question, see e.g. |BLS08[ IKN09] . For more intricate properties, like the regularity of the 
density of states or the analysis of spectral fluctuation boundaries, the difference between 
monotone and non-monotone models is even more striking. 

Nevertheless there has been a continuous effort to bring the understanding of models with 
non-monotone dependence on the randomness to a similar level as the one for monotone 
models. The model of this type to which most attention was devoted so far is the alloy 
type Hamiltonian with single site potentials of changing sign, see e.g. |Klo95t IStoOOt IKlo02t 
IVes02t IHK02t IKV06t IKN09] . More recently also the discrete analog of this model was 
studied in |VeslObt lETVlOt IVeslOat ITVlOt IKrii] . Electromagnetic Schrodinger operators 
with random magnetic field |Uek94l lUekOOl IHK021 IKNNN031 IUek08l IBou09j. as well as 
Laplace-Beltrami operators with random metrics |LPV044 ILPPV08"| ILPPV09"] are other 
examples with non-monotonous parameter dependence. 

A very interesting model with geometric disorder is the random displacement model, 
which exhibits also no obvious monotonicity, cf. e.g. |Klo93j lBLS08t IKLNS] . Another rele- 
vant model (although not defined in terms of a countable family of i.i.d. random variables) 
without obvious monotonicity is a random potential given by a Gaussian stochastic field 
with sign-changing covariance function, cf. |HLMW0T| IUek044 IVes] . 

In this paper we consider a family of Hamiltonians which consists of the Dirichlet Lapla- 
cian on a randomly wiggled waveguide. In this model the dependence of the quadratic 
form on the random variables is neither monotone nor linear. In this respect it is related 
to the random displacement model. Moreover, in our model, the randomness does not enter 
via potential terms, but rather through differential operator terms. A random waveguide 
model has been studied before in |KS00] . There the randomness enters via a variation of 
the width of the waveguide. This type of perturbation leads to a quadratic form which 
depends monotonously on the random variables and is thus structurally different from our 
model. 
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There is a substantial body of literature devoted to the analysis of eigenvalues below 



the essential spectrum of bent, asymptotically straight waveguides, see for instance |ES89 



These eigenvalues are in contrast to the purely absolutely continuous spectrum exhibited 
by an straight waveguide. Thus a local geometric perturbation leads to the emergence of 
discrete eigenvalues. Given this fact, it is interesting to ask whether geometric perturba- 
tions which are ergodic and random lead to dense point spectrum below the continuous 
one, in analogy to the phenomenon encountered for several classes of random Schrodinger 
operators mentioned above. One should point out that in the present paper the local geo- 
metric perturbations of the waveguide are introduced in a somewhat different way than in 
ES89j . 



Let us now describe the main result of this paper. We derive lower bounds on the first 
eigenvalue for a finite segment of a randomly wiggled strip in R 2 . They measure how 
far the eigenvalue may move up, if the vector of random variables moves away from the 
optimal configuration. As an application we obtain a second result: In the terminology 
of the multiscale analysis (MSA) it is a initial length scale estimate for energies near the 
bottom of the spectrum in the weak disorder regime. It corresponds to the induction 
anchor of the MSA. This should be understood as a step towards a localization proof using 
the MSA. If there would be an appropriate Wegner estimate at disposal at low energies 
(which we don't have at the moment) an adaptation of the usual MSA, e.g. as presented 
in [StoOll IGKOll IGK04| would lead to localization. 

Let us say a few words about our methods of proof. It consists of a deterministic 
and probabilistic part. For a finite segment of the waveguide one can use methods from 
asymptotic analysis to estimate the position of the principal eigenvalue of the Laplacian. 
In this situation only a finite number of random variables enters the operator. It is this 
part which requires the weak disorder restriction. The mentioned results can be combined 
with a Dirichlet-Neumann bracketing argument and a large deviations principle to arrive at 
an exponential probabilistic bound on the position of the lowest finite segment eigenvalue. 
Using a Combes-Thomas estimate |CT73l IBCH97t IStoOll IBdMS03] this can be turned into 
an off-diagonal decay estimate on the Green's function, which plays the role of the initial 
length scale estimate in the MSA. 

It is maybe worthwhile to point out some differences to the recent paper |KN09] of Klopp 
and Nakamura which is devoted to the proof of Lifshitz tails for alloy-type Schrodinger 
operators with single site potentials which are allowed to change sign. There are two 
aspect in common between this work and ours: both of them concern the analysis of the 
low lying eigenvalues of finite volume random Hamiltonians and both of them deal with 
non-monotone parameter dependence. There are also two main differences: we are not able 
to give an Lifshitz bound on the integrated density of states for our model, since the global 
disorder coupling constant has to be chosen dependent on the volume scale. If one lets the 
scale tend to infinity the global coupling constant has to go to zero. On the other hand 
we assume no reflection symmetry for the individual perturbations. A crucial assumption 
of |KN09] is that the single site potential obeys such an condition. It allows Klopp and 
Nakamura to use, after some work and ingenious ideas, an effective decoupling between 
different random variables (similarly as in the case of fixed-sign single site potentials). 
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This means that it is only necessary to perform perturbation theory with respect to one 
coupling constant. This aspect of the proof of [KN09] is discussed on page 1134 before the 
statement of hypothesis (H2) there. In our model there is no such symmetry assumption. 
This means that the analysis of the single parameter random Hamiltonian on a unit cell 
with Neumann b.c. does not give us the crucial information which was instrumental in the 
proof strategy of [KN09j . Consequently, we need to analyze the fully interacting model, 
which leads to an eigenvalue perturbation problem with respect to many parameters. 

The paper is organized as follows: in the next section we define rigorously our model, in 
Section [3] we state the probabilistic estimates on the position of the principal eigenvalue of 
a finite segment of a random waveguide and on the exponential off-diagonal decay of the 
associated Green's function, in Section 0] we reduce the proof of the two above statements 
to deterministic bounds on the first eigenvalue, and in the final Section the mentioned 
deterministic estimates are established. 

2. Model 

We consider random quantum waveguides in R 2 , determined by the following data: 
Let (uikjkez be a sequence of independent, identically distributed, non-negative, bounded, 
non-trivial random variables, n > a global coupling constant, / ^ 1 the length of one 
(periodicity) cell of the waveguide, and g E Cq(0,1) a single bump function. The following 
function G:Rx04R determines the shape of the waveguide 

G(xi,u) := y^ j m k g(x 1 - kl). 

Note that kG(xi,u) = G(xi,nu). Together with the global coupling constant k > it 
defines an infinite waveguide as the set 

D KjU : = {x E R 2 | X\ E R, k,G(xi,u) < X2 < kG(xi,u) + n} 

= {x E R 2 | X\ E R, G(xi, kuj) < x 2 < G(xx, ku) + n} = D l KUJ 

Our results are formulated in terms of spectral features of finite segments of the infinite 
waveguide. We define them next. For N E IN and j E Z set 

D K , u (N,j) := {x E R 2 | jl < X\ < (j + N)l, kG(x x ,u)) < x 2 < kG{x 1 ,uj) + n}. 

Denote by T K ^(N, j) the upper and lower part of the boundary of D K)U (N,j), i.e., 

FkA N J) -={x e r2 I J 1 < x i< (j + N)l,x 2 = kG(x u u)} 

U{x E R 2 | jl < X! < (j + N)l, x 2 = kG(xi,u) + vr}. 

The remaining part of the boundary dD K bJ (N, j) \ T K>u (N,j) is denoted by •y KtbJ (N,j). 

Let T-L KjUJ (N,j) denote the negative Laplace operator on D KjUJ (N,j) with Dirichlet bound- 
ary conditions on T KiU (N,j) and Neumann b.c. on 7 KW (iV, j). The lowest eigenvalue of 
HK,u>(N,j) is denoted by X K)W (N,j). If j = we shall use the following shorthand notation: 

D K , U (N) := D^(N, 0), T^(N) := T K>U (N, 0), j K , u (N) := j K:UJ (N, 0), 

U K>U {N) := U K ,UN, 0), \^(N) := X K>U} (N, 0). 
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Similarly as for the infinite waveguide we have D K)W (N,j) = D ltKLU (N, j). Since k > is 
arbitrary we may assume without restricting the model 

(2- 1 ) maxUs-Hqo,/], || ||#'||c[o,z], lb"||o[o,q} = 1. 

Denote the distribution measure of Uk by fi. It will be convenient to think of fi as 
a measure on the semiaxis [0, oo) with support in the unit interval [0,1]. Thus any ujk 
takes values larger than 1 only with zero probability. Then P = /i denotes the product 



measure on the configuration space f2 = x fceZ [0,oo) whose elements we denote by oj = 

Note that by the assumptions on [i the following statements hold for P-almost all u G Q: 
oj G ^°°(Z) and \\uj\\oo '■= su Pfcez \ u k\ = sup fceZ u;fc < oo. On appropriate subspaces of Q we 

define the usual £ p -norms, in particular ||tu||i := ^2 and ||w||2 := ( 

feez fcez 
The randomness of the finite segments D K ^(N,j) arises from only a finite number of 
random variables (u)k)- For this reason it is convenient to have the following notation at 
disposal: For AcZ define the projection map 

Here \a denotes the indicator function of a set A. If A C Z is finite ir A (Q) is contained 
l p {fL) for every p G [1, oo]. We shall use the shorthand notation ||w||a,p := IKa^Hp and in 
the case A = {1, . . ., N} 

IMIiv,p := IMU,p := IKa^IIp- 

3. Probability of low lying eigenvalues and the initial length scale 

estimate 

Here we present estimates on the probability that the lowest eigenvalue of % K) u(N,j) is 
close to one. The first information which one has is that the minimum of the spectrum of 
the Laplacian on an straight waveguide segment %Q )UJ {N,j) is equal to one. This follows 
directly from separation of variables. We shall see later that no operator H K)U) {N,j') has 
spectrum below one. In this sense we can say that one is the minimal spectral value for 
all random configurations. 

To formulate the main result we shall use the following quantities: 

1 f l /j\ ii 3||<7||l 2 ( ,i) 3 llfl'll! 2 (o,«) 

9 9 ~ T / 9{t) dt, c 2 = — , c 3 - 



I J Q aw ' 2P ' 5000/ 7 ' 

Theorem 3.1. Let g and \i as above be given. Let 7 > 34. Then there exists an initial 
scale Ni such that if N ^ Ni the interval 



in 



,c 3 A r ^ 



is non-empty. If N ^ N\ and k G In, then 

(3.1) p(wgO KA n ) ~ 1 < N ~^) < N 1 '^ e~ CiNlh 
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for a constant c 4 > depending only on \x. 

Remark 3.2. The definition of the interval In encodes how our weak-disorder regime de- 
pends on the length scale N. As an instance let us choose 7 = 35. Then it is possible to 
choose k = CsN~^, which means that the allowed disorder regime does not shrink very fast 
for N 00. For larger 7 the behavior is even better. 

Using a Combes-Thomas estimate we arrive at the following estimate on the probability 
that the Green's function, resp. resolvent, exhibits exponential off-diagonal decay for 
energies very close the energy one, i.e. the overall minimum of the spectrum. Note that 
for N large I N C [0,1]. 

Corollary 3.3. Let g, fi, In, 7, N\ and C4 be as in Theorem \ 3.1\ Let k & In C\ [0,1], 
a, (3 ^ 2 and set 

A := {x G R 2 I < xi < a, kG(xi,oj) < x 2 < kG(xx, u) + 71} C D Ki<JJ (N), 
B := {x G R 2 I L - < xi < L, kG(x 1 ,u) < x 2 < kG(xi,u) + n} C D Ki<JJ (N). 
Then we have for any N ^ N\ 

P L G I V A G [1, 1 + 1/(2^)] : - A) _1 Xb|| ^ 2v / iVe(-^ i Sw i )^) 

^ l-iV 1 "^ e" c ^ 1/7 . 

In the typical formulation of an initial scale estimate for the multiscale analysis one 
requires the probability of 

{3 AG [1,1 + 1/(2^)]: \\xa{U k ,u{N)- Xy'xsW > 2VNe~ ± ^W 1 } 

to be bounded by an inverse power of N. Let q G IN and N be such that N {1+q ~ 
exp(A^) for all N ^ jV . Set N 2 := max(Ai, N ). Then for any N ^ N 2 we have 

P G I V A G [1, 1 + l/(2\/iV)] : ||xa(H Kj1 ,(A) - A)^ 1 X b|| < 2ViVe-(-^)) 



>l — N~ 



4. Proof of Theorem 13.11 and Corollary 13.31 



In this section we prove that Theorem 13.11 is implied by the following Theorem and its 
Corollary. We also show how Corollary 13.31 follows. 

Theorem 4.1. Recall that (uikjkez is a sequence of non-negative reals. Let g, I, H KA ,(A), 
A K , W (A) be as in Section^ Recall that 



1 f l 

9 ■= 9 - -j J g(t) dt. 
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«\\ u \\w<^\\MUwjw = w 

Then the estimate 

3 o K 2 \\u\\ 2 Ni2 K 2 \\u\\n,2 



c 2 



AT3 



holds true. 



Corollary 4.2. If k ^ c^N i , then we have for F-almost all u G O: 



r 2 

\ K , u (N)-i> c 2 r 



N 3 ' 

Proof. This follows immediately from Theorem 14. 1[ since ||o;||jv,2 =S; V^V||ca || jv.oo ^ V^V for 
almost all cu. Thus k ^ constant N~^~ implies «||w||2 ^ constant N~ 7 . □ 

The following large deviations principle will be used in the proof of Theorem 13.11 

Lemma 4.3. Let ujf.,k G Z be an i.i.d. sequence of non-trivial, non-negative, bounded 
random variables. Then there exists a constant C4 > depending only on fi such that 



n 



E{u k } 



VnelN : P w - > u k ^ ^ e 



-C4n 



n 
fe=i 



Proof of Theorem \3.1[ For G 2IN and 7 G IN we set Af := K 1 '. Note that we can 
decompose a long waveguide segment into smaller parts. Thus, up to a set of measure zero 
D KiU (N) equals 

U D «AKJ) 

j=0,...,J-l 

1 • 

where J = N/K = K 1 " 1 = N ~~ and |J denotes a disjoint union. According to the 
decomposition of the segment D K}U (N) we introduce new Neumann boundary conditions, 
which decreases the operator. More precisely, we have in the sense of quadratic forms 

j -1 

j=0 

In particular, 

(4.1) \ KtU (N)>min\ KtU (K,j). 

3=0 

The above considerations can be turned into a probabilistic estimate on the position of the 
lowest eigenvalue. Similar ideas have been used e.g. in [MH841 IKSS98] to obtain an initial 
scale estimate. First note that by (14.11) we have the inclusion 

j-i 

{00 G n I A^(iV) - 1 ^ c[j{u6fi \^{K,j) - 1 K-l} 

3=0 
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Since the random variables u^, k G Z are independent and identically distributed, 
^P^weO X K A K J) ~ 1 < ^~ ? ) < A K , w (Ji) - 1 < K~ 

By Corollary 14 . 2 1 and ^ V^IMI^ the following inclusions hold for all k ^ c^K" 1 ^: 

1 f ^2 1 1 Co? 1 1 2 

1 ,-- 3-7/2 1 f llwllifl 1 T i 7 

KJO2 I I A ^VC2 



Denote by Ejcjfc} the expectation value of (any) u k , and choose now k such that 



(4.2) ^ 1 1 i.e. — — s k. 

The upper and the lower bound for k can be reconciled if 7 > 34 and 

2 

2 \ 7-34 

K^K X :- 



E{uj k }c 3 y/C^ 



-2 7 



The last inequality is equivalent to N ^ iVi := (E{wfc}c3- v /c2/2)'i'- 34 . For k satisfying 
we are able to apply the large deviations principle of Lemma 14.31 and thus obtain 

m ( , IMIi 1 r ,i A ™/ , IMIi E{w fc }\ rtK 



K " / \ K 

for K ^ K\. Consequently we have for N ^ Ni 

F (u eQ\ \ K ,„(N) - 1 ^ N 1 -^ e~ C4Nl/ \ 

Proof of Corollary \3.3\ . Since A K)CJ (iV) is the element of ai^K K>u} {N)) closest to one, 
Q' = {weO| \ K ,„(N) - 1 > AT§} = jo, e I dist(l,a(^(A^)) > AT§} 



□ 



wGfi|VAe[l,l + l/(2ViV)] : dist(A, o-(H K , u (N)) > 1/(2VN) 

For w 6 fl' we want to use a Combes-Thomas argument to obtain from the estimate 
on the distance of the relevant energies to the spectrum a decay estimate for the Green's 
function. Such estimates are rather standard for Schrodinger-type operators. However, 
we did not find a specific formulation of this result in the literature which fits exactly our 
situation, so we provide a direct proof, for the convenience of the reader. 

In the following the parameters u G Q', k > 0, and N G IN will be kept fixed. For this 
reason we shall abbreviate in the subsequent calculations 7i K>u (N) simply by %. Define 
the function J: [0,/] -> [0, 00) by setting J(t) = t for t G (1,L-1), J(t) = 3t 2 - 3t 3 + 1 4 for 
t G [0, 1], and choosing the values on the segment [L — 1, L] in such a way that the graph of 
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the function becomes point symmetric w.r.t. (L/2, L/2). Note that J satisfies Neumann 
b.c. at t — and t = L, that it is twice differentiable, and that on the segment t G [0, 1] 

J\t) = 6t- 9t 2 + At 3 ^ 0, J"(t) = 6 - 18* + 12t 2 . 

It thus follows that 

imioc^, Halloo < e, 

and that J is monotonously increasing. For a G (0, 1) we define the multiplication operator 

T a . L 2 (D K>UJ (N)) L 2 (D K>UJ (N)), (T a f)(x) : = e aJ ^f(x) 
and another operator 

7/ d ,2 V 



V a := -2a J' 



dxi 



a 2 J' - a J" 



which will turn out to be an 'effective perturbation'. A direct calculation shows that 
T-*% = I, T-aUTa = U + V a , T- a {H - X^Ta = (H + V a - X)~\ 

provided that A is in the intersection of the resolvent sets of the two operators % and 
T-L + Va- We shall identify a range of values for a such that the last condition holds and we 
get even an explicit bound on the norm of {% + V a — A)" 1 . In these considerations we shall 
keep in mind that A is close to, but larger that one, in particular A G [1,2]. We denote by 
5 = dist(cr("H), A) the distance of the spectrum of "H to A. Let us first estimate 

d 



(4.3) 



(«ll^llL + ll^lloo)||(^-A)- 1 ||+2||J'| 



dxi 



(n - a)- 1 



< a 



ffa + 6 



A 



Here we have employed that 







dx\ 



(H-X) 



-i 



A 

6 2 



which follows directly from the obvious relations 



ll Vn llL(D^) - A IMli 2 (D K ,.) = tf,u)l*{D K , u ), 

where u:=(H- A)" 1 /, / e L 2 (D K>UJ ). 

The right hand side in H4.3[) is bounded by 



L 2 {D K 



L 2 {D K 



121 5 
165 + 2 



A + 5 



a 

*12- 



since 6, a G [0, 1] and A G [1,2]. Now choose a = £ \±. Then \\V a (H - A)" 1 )] «C 1/2 
and thus the norm of the Neumann series satisfies 

\\(n-xrn _ 2 

l-i <T 



\\{H + V a -X)- l \\ ^ 
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Now choose a, (3 ^ 2 and set 

A := {x G R 2 | ^ xi ^ a, kG(x x ,uj) < x 2 < kG(x x ,u) + n} C D KjU (N), 

B := {x G R 2 | L - f3 < xi < L, kG(x x , oj) < x 2 < kG(x x ,u) + 7r} C D KiUJ (N). 

For any normalised vectors <j),ip £ L,2(D KjUJ (N)) we have 

KivlxA.r-aCw-A)- 1 ?;^^!)! ^ ||7i„(«- a)- 1 ^! = iKft+Pa-Ani < 

Due to the monotonicity of J and the positivity of the integral kernel of (% — A) -1 we are 
able to estimate 

|(MXA, 71.(7* - A) _1 7^Xb|0|)| ^ e a ( J ( L -«" J ( a » |(to, - A)"W>| • 

Note that by the choice of the values of a, /3 and the function J we have J(L — 0) — J (a) = 
L — (3 — a = dist(A, B). Bringing the exponential term on the other side we obtain 



| (to, (ft- A) ~W} | < 



2 _odlst(A,B) _ g dfat(A,B)< 

5 ~ 5 



24 



Now fixu£ff = {wGO A« )W (iV) - 1 > jV-a J and A G [1, 1 + ^-]. Then 5 ^ ^= and 
thus 

■ . i , i 2 dist(A,B)« , dist(A.B) 

| (to, (7£ - A)~ < ~± — < 2ViVe «v* . 

Since we have by the estimate (13.11) the bound P(fl') ^ 1 — A^ 1- ^ g-^N 1 / 7 f or N ^ N x , we 
conclude that 

(4.4) p(wen|VAe[l,l+ 1/(2VF)] : ||xA(ft«^(iV) - A)- 1 ^!) < 2v / iVe"^W# 1 ) 

□ 

5. Deterministic lower bounds 

In this section only finite segments of the waveguide will be relevant to us. Recall that 
due to the fact that the support of the measure \i is contained in the interval [0, 1], for 
almost all u G Q the bound \\u>\\oo ^ 1 holds. In this section an ^-normalisation will be 
better suited, and in fact possible since only finite waveguide segments, and thus only finite 
families of random variables u x , . . . are involved. 

We fix iV G IN and set for all j G {1, . . . , N} 



Pj := KUj,9j := Pj/e, where e := ( ^~]p 



N \ 2 

2 



Observe that the vector 9 = {9i\f =1 is normalized in the sense 
(5.1) 9\ + ... + 9 2 N = l. 
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Let / ^ 1, g : R ->• E, g E R, G : R x R be as in Sections [2] and H We set L := Nl 
and consider in the following the restriction of G to the set [0, L] x ( x jLi [0, oo)). The 
restricted function will be again denoted by G, and the finite vector (ooi, . . . ,oon) will be 
denoted by oo, by slight abuse of notation. In terms of the new parameters we have 

JV 

kG(x±,uo) = G(xi, koo) = G(x±, p) = G(xi, e9) = eG(xi, 9) = e 9jg(x\ — jl). 

For the subset of R 2 forming the finite segment of the waveguide we shall use in this section 
the notation 

Il p := {x G R 2 | < xi < L, eG(xi, 9) < x 2 < eG(xi, 9) + 71} 

Let us point out that U p = -D||p|| 00l p/||p|| 00 (-/V') in the notation of Section [2] and that n = 
{xGR 2 |0<xi<I,0<i 2 < tt}- Similarly as before 

T p := {x e R 2 I x 1 e (0, L),x 2 — eG{x u 9)} U {x : xi e (0, L),x 2 = eG(xx, 9) + n}. 

indicates the upper and lower parts of the boundary of U p and the remainder of the 
boundary dH p \ T p is denoted by 7 p . By *H p we denote the negative Laplacian in U p with 
Dirichlet boundary conditions on T p and Neumann ones on 7 p and by A(p) the lowest 
eigenvalue of H p . 

In what follows (•, •) and || • || denote the scalar product in L 2 (IIo) and the associated 
norm. We repeat here the statement of Theorem 14.11 in the new parametrization. 

Theorem 5.1. For 

3 1 

(5-2) e< 5000 Nli ^^ 

the estimate 

3 e 2 

holds true. 

The rest of the section is devoted to the proof of the theorem. 



First we transform the Laplacian on the wiggled strip segments into a differential opera- 
tor with variable coefficients on a rectangle. It is then possible to treat the latter operator 
as a perturbation of the pure Laplacian on the rectangle. We introduce the coordinates 
£ := £2), £1 : = { 2 := x 2 —eG(xt). The mapping u(x) t-s- w(£i,£2 +eCr(£i)) is a uni- 
tary operator from L 2 (U p ) to L2(IIo). Let ip p be the normalized eigenfunction associated 
with \(p). It is the unique solution of the boundary value problem 

(5.3) (-A € - eQ p )i> p = X(p)^ p , £ E n , 

(5.4) ^o = 0, xeT , 7^ = 0, ee 7 o, 

»£i 
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Q P ■= -2G'- 



Hereafter the prime denotes the derivative w.r.t. £i- 

For the perturbation theoretic estimates we are aiming for we need some control of the 
resolvent set. The first two eigenvalues of "Ho are 1 and 1 + tt 2 /L 2 . The eigenfunction 
corresponding to the eigenvalue one is 



^ :n -^n , ^o(0 := y sin£ 2 . 
It turns out that it is more convenient to work with the modified resolvent 



fto(A) := (Ho - A)" 1 - 



which is a bounded self-adjoint operator for any 

(5.5) A G B„ 2/(2L2) (1) c C. 

The operator-valued function B^/(2L 2 )(X) 9 A 4 7^ (A) is holomorphic, cf. |Kat66[ Ch. 
V, Sec. 3.5] 

For an arbitrary / G L 2 (n ) set f '■= f — (f, i>o)i>o- Then the function u := TZo(X)f = 
(Ho — A) -1 / solves the equation 

(Ho - X)u = (Ho - A) (Ho - A)- 1 / + (Ho - A)^^Vo = /, 



(5.6) 

and is orthogonal to ipo m ^(U-o)- 



A- 1 



Lemma 5.2. For f G L 2 (n ), u := 1Z (\)f and A G B^/^l 2 )^) C C we haveu G iy 2 2 (Il ) 
and 



(5.7) 



2L 2 


u\\ ^ - 




TV 2 


du 













IVttll ^ 



7L 2 



7T 



25L 2 



7T Z 



<9 2 W 



AIL 2 



7T 



Proof. The vector u = (Ho — A)" 1 / is in the range of ("Ho — A)" 1 and thus in the Sobolev 
space W / 2 2 (n ). The first inequality follows from |Kat66[ Ch. V] and the fact that our choice 
of A is separated by at least n 2 /(1L 2 ) from the next spectral value. 
Let us prove the second estimate. We begin with the obvious identity 



(5i 



\(f^o) 2 \ ^ 
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We multiply f 15 . 6 j) with u and obtain || Vw|| 2 = A||w|| 2 + (/, u). This identity, (15. 5p and the 
first estimate in (15 .7D yield 



|Vu|| 2 ^ lAllld 



(5.9) 



at^ or 2 
lnll^^|A|||/|| 2 + ^||/|| 2 



7T J 



7P 



4L 4 



71 1 



1 + 



7T 



7T 



2L 2 2L 2 



4L 4 



71 J 



(1 + ^) 



that proves the desired estimate. 

Since u e W^IIo), the function v := is a generalized solution to the boundary value 
problem 

df 

-A^ = A^ + -f, een , v = o, £edii 

which is obtained by differentiating equation (15. 6p . We multiply the last equation by v, 
integrate by parts, and obtain 



X\\v\ 



We employ (I5.5p . (I5.8p . and Young's inequality to estimate 



|Vd| 2 ^ 1 



7T 



2L 2 

and hence by (15. 9 p 



, du 



Wv\\ 2 < 2 



dv 



< 1 



7T 

2L 2 



|Vn| 



iVfl 



7T 
L 2 



|Vn|| 2 + 



(5.10) 



Ti- 



ll 



^ + 1)^2 + ^+-^- 

^(4( 7 r 2 + l)(2 + 7 r 2 ) + 7 r 4 )^||/|| 2 ^^ 

This implies the penultimate inequality in (15. 7p . 
We rewrite (15 .6p as 

d 2 u d 2 u 

and see that the estimates (15. 5p . (I5.8p . (I5.10p . and the first estimate in (15 .7p yield 



d 2 u 



d 2 u 



+ |A| Hull + 



^ v / 614 + 2 1 



Ti- 



ll? 

^V / 614 + 2(l + 7T 2 ) +7T 



7T 



L 2 



Ti- 



7T Z 



47L 5 



7T 



which proves the last claim in (15. 7p . 



□ 
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The inequalities in (I5.2p can be understood as bounds on the norms of certain operators 
products: 



(5.11) 



l^o(A)|| < 
v|^ (A) 



2L 2 



7T- 



25L 2 



7T- 



|V^ (A)|| ^ 
d 2 



7L 2 



7T 



47L 2 



By the normalization (12. ip we have 



(5.12) || 

9\\l 2 (o,i) — IMIl 2 (cm) 



T x ( J g(t)dt)\\\g\ 



7T Z 



L 2 (0,0 



^ I. 



This implies together with inequality (15. 2p 

3/2 3 1 3 

^ ' ' £ ^ 5000L 7 ^ 5000 L 5 N 2 ^ 5000' 

We use (15. ip and (12 .ip and Cauchy-Schwartz inequality to establish 

\\Q P u\\ ^2 

for any u G W^Ho). For A satisfying (I5.5p . and e as in (I5.13p . Lemma [5T21 and the last 
bound imply 



d 2 u 




d 2 u 




du 




+ 


da 


+ 


96 



(5.14) 



1 0AT 2 

\Q p K (\)\\ ^ e\\Q p Ko(\)\\ < 



TT* 1257T 2 

Here the norm is understood as the norm of an operator in L 2 (n ). 
Let us show that A(p) satisfies (15.51) . Indeed, 



< 1. 



(5.15) 



difjp 



eG 



96 



2 

+ 


dip p 


2 


dip p 




96 


96 



^ 1. 



By the minimax principle and (15. ip . (12. ip the eigenvalue X(p) can be estimated from above 
as follows 



(5.16) 



di/j 



eG 



,9^o 



96 96 
Two last estimates and (I5.13P imply (15. 5p for A(p). 



2 


9^o 


+ 


96 



c 1 



Next we shall do some perturbation theory for linear operators to be able to identify the 
dominating contributions. Using the proved fact, (I5.14p . and proceeding completely as in 
|Gad02j . |Bor06t Sec. 4], one can show that for the considered values of e the eigenvalue 
A(p) solves the equation 

(5.17) A(p) - 1 = - £ <(I-eQ p ^o(A(p)))^ 1 Q p ^ ,^o>, 

where I is the identity mapping. A direct calculation shows 

(I - sQpK Q (\(p)))- 1 = I + sQ p n (\(p)) + e\\ - eQpK^p)))- 1 (Q p K (\(p))) 2 . 
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We substitute this relation into (I5.17p . 

Hp) - 1 = - e(Q P ?Po, ^o) - e 2 (Q p n (\(p))Q P ij , ^o) 
5 ' 18 -^<(I-^Qp^o(A(p)))^(Q p ^o(A(p))) 2 Q^ ,^o>. 
We first consider the coefficient of e 2 . For the sake of brevity set 

Integrating by parts, we obtain 

= " 2 \ G V ^/ + \ G w ' WJ + £ \ {G )V 'W/ 

dip \ 



G / X(A(p))Q p ^o,^>-^((G") 2 ^o(A(p))Q p ^o,^o) 



which equals 

G"n (X(p))G"^, ^) + e{G"TZ {\{p)){G'f^ ^) - e((G') 2 n (\(p))Q p ^ 

Since 7^ (A) is a (special) rank one perturbation of (Ho — A) -1 , cf. (15. 6p . it inherits the 
resolvent identity 

Ko(\) - 72o(a0 = (A - h)TZq{p)TZq{\). 
Now we substitute two last identities into ( 15.18p . 

A(p) - 1 =-e(Q p ^^ )-e 2 (G"n (l)G"^^) 

+ £ 3 ((G / ) 2 ^o(A(p))Q^ ,^o) 

(5.19) -e 3 (G'X(A(p))(C7') 2 ^o,^} 



9C: 

-£ 2 (A(p)-l)(C7^ (l)^o(A(p))G^ 



-£ 3 <(I-£Q^o(A(p)))^(Q^ (A(p))) 2 Q^o,^o>. 

It will turn out that the two first terms on the right hand side of this identity are positive 
and dominate all the remaining terms. The first term is of the order e 2 /L. However, a 
cancelation occurs such that the leading contribution turns out to be of order e 2 /L 3 . A 
direct calculation shows that 

(5.20) (Qp^Po^o) = -jh'Wlm- 
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It is more complicated to evaluate the second term. We denote 



w := TZ (1)G 
Since this function solves the equation 

(H — l)w — G 
it is natural to consider the auxiliary function 
(5.21) v:=w + G 



96' 
dijj 



96 

Using the normalization (15. ip . separation of variables and 



d£ 2 llL 2(lVr) 



1/L, 



we get 
(5.22) 



n"T> (-\\n" d ^ d ^°\ In" d ^°\ 1 n 'i|2 , / ^" d ^°\ 

G n ° (1)G W WJ = \ Gw 'WJ = i h L >™ + v> G WJ- 



96 



96 



The first term on the right side leads to a cancelation of (15.201) . It is clear that the function 
v solves the boundary value problem 

-Av = v, £ e n , ^- = o, £ g 7o, 

96 



J-, 6e(0,L), 6 = 0, 



J-, 6e(0,L), 6 = vr. 



We multiply the equation in this problem by G^^ and integrate twice by parts. Due to the 
separation of variables a direct calculation shows that part of the boundary contributions 
vanish and the equation simplifies to 



v,G 



96 




7rL 



-<v,G 



G(6 

>,dip 
96 



dv 
96 

v,G 



dv 

6 =7r 96 
dip 



d6 



96 



Thus 



(5.23) 



v,G 



96 



dv 

^ 2 =7T 96 



+ 



6=o 



d6- 
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We shall now expand the functions in terms of the eigenvectors of the transversal Laplacian 
with Dirichlet b.c. and the longitudinal Laplacian with Neumann b.c. It follows from the 
equation (H - \)w = G"g| that 

[~2~ n 2 mG n (p) cosnnL' 1 ^ sin2m£ 2 
- 2^ 4m 2 -1 (4m 2 - 1)L 2 + ix 2 n 2 ' 

n,m=l v ' 

where the coefficients G n (p) are defined by the identity 

oo 

G(£i,9) = Y,G n (p) cosvrnL- 1 ^. 

n=0 

We substitute the expansion obtained and (I5.2ip into (I5.23p . 



9£ 2 / ' 4m 2 — 1 (4m 2 - 1)L 2 + 7r 2 n 2 

2/^2 , 



32^ n 2 G 2 n (p) 32 ^ 2 64 2 

> 3 ^ 3L 2 + tt% 2 > 3(3 + tt 2 )L 2 ^ n[P) > (9 + 3tt 2 )L 2 

n=l v n=l v 

Here we have also employed that due to ( 15.121) and the normalization (I5.20p 

ELG^p) _ .. ||2 _ ||^i|2 j n 2 

— - — - ||Lt - Ct || L2 ( 0iL ) - II^IIlzCo.l) ~ LLt o 



N 

|2 



9\\U,i)- T \^i) 9(t)dt 




> lblli 2 (o,o - y ( / 9(t)dt ] = \\g\\l 2m . 
Now it follows from (OP]) . (15321) that 
(5.24) - e(Q p ^ - , 2 (g'X(1)G''§|, §|) > |^||<7|| W 

Using (EE}, f lBTTTT) . flBTTBl . (EH, llf^ll = 1, and the inequality 

HQ^oll «S l + £^2, 
we estimate the remaining terms in the right hand side of (I5.19P by noting that 

e 3 |((C7 / ) 2 ^o(A(p))Q^o,^o)| < ^ 



Ae 3 L 2 4e 3 L 4 



£ 3 



G"K {\{p)){G')% 



A 2 / <^ 



96 



2e 3 L 2 2e 3 L 4 
^ — ^- ^ 



7T 2 7T 2 
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To estimate the penultimate term (15.19P observe that by (15.15p . (15 . 1 6[) and (15 . 1 3[) 



A(P) - 1) 



35 



3s 



5000L 5 iV 2 5000L 5 
since N ^ 1. The last estimate, (I2.ip . and the first inequality in ( 15. lip yield 



(X(p)-l)(G"7Z (m (X(p))G 



3e 3 /2L 2 \2 



5000L 5 V vr 2 



dtpo 



L 2 (n ) 



3e 3 3s 3 L 4 



1250vr 4 L 1250vr 4 ' 



We use (15.141) and the estimate 



to arrive at 



s 3 \((l-sQ p n (\(p))r\Q p n (\(p))) 2 Q^ ^ )\ < 250 • ^ L l e 3 . 



125tt 4 - 8tt 2 



In view of ( 15. 19ft . ( 15. 24ft it leads us to the estimate 

64£ 2 



(5.25) 
We note that 



Hp) - 1 > 



(9 + 3vr 2 )L 3 



||^|| 2 2(0i0 -225LV. 



64 



225 3 
9 + 3tt 2 ~ 5000 > 2 

and use assumption ( 15. 2ft to bound ( 15. 25ft from below by 



64e 2 



(9 + 3vr 2 )L 3 
which completes the proof. 



\9\\l 2 (o,i) 



- 225LV 



3||?| 



Li(o,i) . 3 , 2 e" 



5000L 7 



2 iiy 11^(0,0^3' 
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